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Abstract. In this paper we consider query versions of visibility testing and visibility 
counting. Let be a set of n disjoint line segments in M? and let s be an element of S. 
Visibility testing is to preprocess S so that we can quickly determine if s is visible from a 
query point q. Visibility counting involves preprocessing S so that one can quickly estimate 
the number of segments in S visible from a query point q. 

We present several data structures for the two query problems. The structures build 
upon a result by O'Rourke and Suri (1984) who showed that the subset, Vs{s), of M? that 
is weakly visible from a segment s can be represented as the union of a set, Cs{s), of O(n^) 
triangles, even though the complexity of Vs{s) can be Q{n'^). We define a variant of their 
covering, give efficient output-sensitive algorithms for computing it, and prove additional 
properties needed to obtain approximation bounds. Some of our bounds rely on a new 
combinatorial result that relates the number of segments of S visible from a point p to the 
number of triangles in Use5 ^sis) that contain p. 



1 Introduction 

Let S" be a set of n closed line segments whose interiors are pairwise disjoint. Two points 
p,q G ]S? are (mutually) visible with respect to S if the open line segment pq does not 
intersect any element of S. A segment s G S is visible (with respect to S) from a point p if 
there exists a point q £ s such that p and q are visible. If two objects (points, segments) A 
and B are visible (with respect to S), then we say that A and B see each other (w.r.t. S). 
In this paper we consider the following two problems: 

Problem 1 (Visibility testing). Given a query point p and a segment s £ S, determine if 
p sees s. 

Problem 2 (Visibility counting). Given a query point p, report the number of segments of 
S visible from p. 

For a point p G M^, the visibility region or visibility polygon of p (w.r.t. S) is defined 
as (see Figure [T}a): 





^s{p) = {(/ G : p and q are visible (w.r.t. S)} . 
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The visibility region of a point is star-shaped, has p in its kernel, and has size 0{n). It can 
be computed in 0(n log n) time by sorting the endpoints of segments in S radially around 
p and then processing these in order using a binary search tree that orders segments by 
the order of their intersections with a ray emanating from p [H [18] . (Equivalently, one can 
compute the lower-envelope of S in the polar coordinate system whose origin is p.) Because 
Vs{p) is star-shaped with p in its kernel it is easy to determine if a query point q is contained 
in Vs{p) in O(logn) time using binary search. In this way, one can consider Vs{p) as an 
0(n) sized data structure that can test, in O(logn) time, if a query point q sees p. 

For a segment s G S", the visibility region of s (with respect to S) 

Vs{s) = V5'(g) = {p G : s and p are visible (w.r.t. S)} 

q£s 

is the set of points in that see (at least some of) s, see Figure [T}b. Unlike the visibility 
region of a point, the visibility region of a segment is a complicated structure. For a segment 
s, Vs{s) can have combinatorial complexity il(n'^) and \^(s) can have Q{n^) connected 
components [HI Figure 8.13] [101 Lemma 12], see also Figure [2] 

More troublesome than the worst-case complexity of Vs{s) is that there exist sets S 
of n line segments where, for most of the elements s £ S, the complexity of Vs{s) is r2(n^). 
Therefore, explicitly computing Vs{s) and preprocessing it for point location does not yield 
a particularly space-efficient data structure for testing if a query point p sees s, even if s is 
a "typical" (as opposed to worst-case) element of S. 

In this paper we propose efficient data structures that use an old result of Suri and 
O'Rourke [18] which shows that Vs{s) can be represented as a set of 0{n?) triangles whose 
union is Vs{s). We define a variant of their covering, give efficient algorithms for computing 
it, and prove additional properties of the covering. In particular, we define a covering Cs{s) 
of Vs{s) by triangles. We prove that for a randomly chosen s £ S, the expected size of 
Cs(s) is 0{n). This, of course, implies that | Use5 ~ O(n^). Additionally, if we 

define C{S) = \Js&s ^s{s), then we prove that the number of triangles of C{S) containing 
any point p is a 2-approximation to the number of segments of S visible from p. 

Applications of these results include efficient data structures for testing if a query 
point is contained in Vs{s) as well as efficient data structures for estimating the number of 
points of S visible from a query point. In order to express our results more precisely, we 
need some further definitions. 

1.1 Visibility Graphs and Extended Visibility Graphs 

The visibility graph VG{S) is a graph whose vertices are the 2n endpoints of the segments 
in S and in which the edge pq exists if and only if the open line segment with endpoints 
p and q does not intersect any (closed) segment in S. (see Figure [3ja). It is well-known 
that the number of edges m of VG{S) is in 0{n^). Ghosh and Mount [12] give an optimal 
0(n log n -|- m) time algorithm to compute the visibility graph of a set of n disjoint line 
segments. Here, and throughout the remainder of the paper, m = m(S) is the number of 
edges of VG{S). 
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Figure 1: (a) The visibility region for a point and (b) The visibihty region of a Une segment. 



s 

Figure 2: An example of a set S where ^^5(5) has complexity O(n^). The 0(n) segments in 
the center define visibility graph edges whose extensions intersect in Q{n'^) points. 



Assume, w.l.o.g., that no segment in S is vertical, so we can say that a point p is 
above a segment s G S if p is above the line that contains s. Assume, furthermore, that 
S contains four segments that define a rectangle that contains all the elements of S in 
its interior. The first assumption can be ensured by performing a symbolic rotation of S. 
The second assumption is only used to ensure that all visibility regions that we discuss are 
bounded. 

The extended visibility graph EVG{S) is obtained by adding 2m edges and at most 
2m vertices to VG{S) as follows (see Figure [3jb): For each (directed) edge uv in VG{S), 
extend a segment Cuv from v in the direction until it intersects an element of S at some 
point w. If not already present, then add the vertex w to EVG{S) and add the edge vw to 
EVG{S). The extended visibility graph can be computed in 0(n log n + m) time using the 
visibility graph algorithm by Ghosh and Mount [121 . 

The union of the edges of EVG{S) and the segments in S form a 1-dimensional set 
whose removal disconnects into a set of 2-dimensional regions. This set of 2-d regions is 
known as the visibility space partition, VSP{S) of S. The regions of VSP{S) are important 
because for any region R G VSP{S) and for any p,q £ R the set of segments of S visible 
from p is equal to the set of segments of S visible from q. The region of VSP{S) that 
contains p determines all the combinatorial information about Vs{p)- 

Note that VSP{S) is defined by ©(n^) fines, rays, and segments and tfierefore has 
worst-case complexity O(n^). 
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Figure 3: The visibility graph and the extended visibihty graph of a set of seven hne 
segments. (Segments are bold, graph edges are dashed.) 



1.2 Previous Work 

There is a plethora of work on visibility in the plane. This section discusses only some of 
the work most relevant to the current paper. 

The visibility space partition is bounded by a subset of the O(n^) lines induced by 
pairs of endpoints in S. The VSP{S) has complexity 0{m?) where m is the number of edges 
in VG{S) and can be computed in O(m^) time after constructing VG{S) using standard 
algorithms. 

By preprocessing VSP{S) with a point location structure and augmenting the re- 
gions of S with appropriate information, one obtains an 0{m?') size data structure that can 
answer visibility testing queries and visibility counting queries in O(logn) time. 

If the segments of S are the edges of a simple polygon then Bose et al. ^ and 
Guibas et al. [13j show that the complexity of VSP{S) is only O(n^). In this case, this 
immediately solves the two problems using a structure of size O(n^). Aronov et al. [3] give 
a data structure that reduces the space to O(n^) but increases the 0(log n) query time term 
to O(log^n), again for the case where segments of 5 are the edges of a simple polygon. 

Pocchiola and Vegter [T^ give an 0(m) space data structure, the visibility complex, 
that can compute the visibility polygon Vs{p) from any query point p in 0(mp log n) time, 
where nip is the complexity of Vs{p). When the segments of S define a polygon with h 
holes then Zarei and Ghodsi [19] give an O(n^) space data structure that can compute 
the visibility polygon Vs{p) in 0(mp log n) time and the query time of their structure is 
0(min{/i, TTip} logn + mp), which improves the query time of Pocchiola and Vegter when 
h <^ nip. 

Motivated by the computer graphics problem of estimating a priori the savings to 
be had by applying a visibility culling algorithm, Fischer et al. |10| [TT] give approximation 
algorithms for Problem 2. They present two approximation data structures for visibility 
counting. One structure uses a (r/m)-cutting [151 Section 4.5] of the EVG{S) to obtain a 
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data structure of size 0{{m/r)'^) that answers queries in O(logn) time and approximates 
the visibihty count up to an absolute error of r. Another structure uses random sampHng 
to obtain a data structure of size {m? log*^^"*^^ t^)/^-, that has query time £log*^^^^ n, and that 
approximates the visibihty count up to an absolute error of 5n for any constant 5 > 0. 
(Note that 5 affects the leading constants of both the query time and space requirements.) 



1.3 New Results 

In the current paper we revisit O'Rourke and Suri's proof that, for any s € S, there exists 
a set Cs{s) of 0{ms) triangles whose union is Vs{s), where is the number of edges 
of EVG{S) incident on s. We show that this covering has the additional property that if 
we take the 0{m) size set C{S) = UseS'C'5(s) of triangles, then the number of triangles 
containing any point p G is a 2-approximation to the number of segments of s that are 
visible from pj^ 

These triangle-covering results have several applications that are obtained by storing 
the resulting triangles in a layered partition tree. Here, and throughout the remainder of 
the paper, e > is a constant that can be made arbitrarily small. To reduce clutter, we use 
the notation Oe{f{n)) = 0{f{n)'rf). 



1.3.1 Visibility testing 

By storing the elements of Cs{s) in a partition tree, we obtain, for any k with < k < m^, 
an 0{k) space data structure that can test, in Oe{ms/Vk) time, if a query point p is 
contained in Vs{s). Barring a major breakthrough on Hopcroft's Problem [8], this result is 



likely only a factor of 0{n'') from the optimal. See Section 3.1 

For comparison, the best previously described structure for this problem, as used 
within the results of Fischer et al. [lOllllj . has size 0{m1/tj and answers queries in O(^logn) 
time, where ^ > 1 is a space/time tradeoff parameter of the data structure. Taking ^ = ^/n 
yields a space of 0{ms ) and a query time of O(y^mjlogn). On the other hand, taking 
k = rri'J'^ in our data structure yields an OiraJ'^) space data structure with query time 
Oe(my'). 



1.3.2 Visibility Counting — Relative Approximation 

By putting all the triangles of C(S') into a partition tree, we obtain a data structure that 
can 2-approximate the number of segments of S visible from any query point. For any k 
with m < k < m?, this structure has size 0{k) and answers queries in time 0{m/ ^/k). The 
structure returns a visibility count m! that satisfies rup < m! < 2mp. See Section 



3.2 



'^In fact, O'Rourke and Suri's covering is a 3-approximation. The slightly modified version we describe 
in this paper is a 2-approximation. 
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1.3.3 Visibility Counting — 



Absolute Approximation 



Using a selective random sampling of the segments in S, we obtain a data structure of size 
Oe{{c'm/n){cn)'^) = Oein^^") that approximates the number of segments of 5 visible from 
any query point in time Oe(c(w./n)(^/^)(^~°)) = Of:{cn^^^'^^^^~°'^), for any given constants 
c, (5 > and < a < 1. With probability at least 1 — n^i^'^'^'^/"^p) ^ the structure returns a 
value m'p such that nip — c/n — 6n < m'p < rup + 6n. This data structure is described in 
Section 13. 3.11 



Using random sampling in a different manner, we obtain a space versus query time 
tradeoff. For any k with m/n < k < (m/n)'^, we obtain a structure of size 0^{k) and query 
time Oe{m/{ny/k)). This structure returns a visibility count m'' that satisfies nip — 6n < 



m!' < 2m„ + 6n. The details can be found in Section 3.3.2 



(a) 




(c) 



Figure 4: The algorithm for covering Vg{s) with triangles processes the events at (a) pi, 
(b) p2 and (c) p^. Active edges are shown in orange and triangles in the covering are shown 
at the time they are added to the covering. 

These results compare favourably with those of Fischer et al. . Their cutting- 

based data structure, with parameter r = 5n, gives an absolute error of 5n, uses space 
0((m/n)^) and has a query time of O(logn). Their random sampling-based data structure, 
with parameter i = -y/n, gives a data structure of size (m^ log*^*-"*^^ n)/-^/n with query time 
y^log'^^^^ n. 

The remainder of the paper is organized as follows: Section [2] proves results on 
covering visibility regions with triangles. Section [3] applies these results to obtain new 
results on visibility testing and counting. Section |4] summarizes and concludes with open 
problems. 



2 Covering Vs{s) 

In this section we give an algorithm for covering the visibility region Vs{s) with a set 
Cs{s) of triangles. The resulting covering is similar to the covering given by Suri and 
O'Rourke [18], the main difference being around the triangles adjacent to the endpoints 
of s. However, our exposition, and our algorithm for computing C{S) are more s-centric. 
This leads to efhcient output-sensitive algorithms for constructing Cs{s), rather than the 
worst-case optimal O(n^) algorithm obtained by Suri and O'Rourke [18]. The number of 
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triangles used in Cs{s) is bounded by 0{ms) where nis is the number of edges of EVG{S) 
that are incident to s. 

We will show how to cover the portion V^(s) C Vs{s) in the halfplane bounded 
from below by the supporting line of s with a set Cg{s) of triangles. The complementary 
part Vg {s) = Vs{s) \ Vg'{s) can be covered with a set Cg{s) using a symmetric algorithm. 

The covering algorithm works by sweeping a point p from left to right along the 
segment s. Events in this sweep occur at the vertices pi, ■ ■ ■ ,Pm'^ of VSP{S) incident on s, 
in their left to right order, so that pi and Pm'^ are the left and right endpoints, respectively, 
of s. 

Let e be an edge of V^{p) that is collinear with p and such that the interior of 
Vg'{p) is to the right of e. We call such an edge an active edge of Vg~{p). Active edges are 
important because, as p moves to the right, they uncover regions of M? which may not have 
been previously visible. See Figure |4ja. 

Let q be the lower endpoint of an active edge e and note that q is an endpoint of 
some segment in S. Consider what happens to e as the viewpoint p moves left to right along 
s, but does not cross any edge of EVG{S) collinear with q. As p moves left to right, the 
edge e remains collinear with p and q and sweeps over a triangle Ag whose lowest vertex is 
q. This continues until the point p reaches an edge of VSP{S) incident on q. See Figure |4|b. 

Algorithmically, the cover C^(s) is constructed as follows: Initially p = pi'is the left 
endpoint of S. We compute the visibility polygon Vg{p), whose boundary is a sequence of 
2mp = 0{n) edges that alternate between subsegments of the elements of S and segments 
collinear with p and an endpoint of an element of S. This polygon can be covered in 
a natural way with nip non-overlapping triangles, each of which has p as a vertex (see 
Figure |4}a). These iJip triangles are added to C^{s). After computing Vg{p) we identify 
its active edges, and with each active edge e we store the value start(e) = pi. 

Next, we sweep p from left to right, pausing at the vertices P2-, ■ ■ ■ ,Pm'^ as we go. 
Upon reaching a vertex pi, we process the edges of EVG{S) incident on pi one at a time. 

Let e' be an edge of EVG{S) incident on pi. If e' is collinear with an active edge e of 
Vg{p) then we generate a new triangle Ae for Cg{s). The lowest vertex of Ae is the lower 
endpoint q of e. Ae is bounded by two lines £±,£2, both of which contain q, and where ii 
contains pi and £2 contains start (e). The third side of Ag is bounded by the segment in S 
incident on e and furthest from pi. See Figure [4]b. 

Finally, the visibility polygon Vg{p) is updated in the neighbourhood of e, which 
possibly creates up to two new active edges incident to q. Each new active edge / is marked 
as active and we set start(/) = pi. The exact nature of this update depends on the relative 
locations of the two segments that define e'. The three possible cases are illustrated in 
Figure [7} 

Note that an important event, but which requires no special handling, occurs at 
the right endpoint of s when p = Pms- In this case, each active edge of Vg{p) generates a 
triangle that is added to the set Cg{s). See Figure |4jc. 

^For segments in sufficiently general position, pi, 1 < i < nis will be incident to only one edge of EVG{S), 
but the covering algorithm does not require this. 
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• • • « • » 

start(e) p' start(e) p' 

(a) (b) 



Figure 5: Proving that C^(s) covers (s). 

We now prove the correctness, construction time and approximation bound of the 
above algorithm. 

Lemma 1. Let Cg{s) be the set of triangles generated by the above algorithm. Then 
UC^(s) = Vg'{s) and \Cg{s)\ < mg where nig is the number of edges of VSP{S) inci- 
dent on s. 

Proof. To prove the bound on the size, first observe that the initial visibility polygon V^{pi) 
has size that is bounded by the degree of pi in VSP{S). Furthermore, at each event point 
Pi, i > 1, the number of triangles added to Cg{s) is at most the number of edges of VSP{S) 
incident to pi. Therefore, the total number of triangles in Cg{s) is at most the number of 
edges of VSP{S) incident on s. 

The fact that UCg{s) C follows immediately from the easily verifiable fact 

that each triangle added to C^{s) contains only points visible from some point on p £ s. 
In particular, for any point r in the triangle Ag that is added to Vg~{s) when processing pi, 
there is a point q in the subsegment of s between start (e) and pi that sees r. 

To prove that Cg{s) covers V^{s), consider a point r G V^{s). If r is visible 
from pi then r is contained in one of the triangles added during the initialization of the 
algorithm. Otherwise, there exists some point p' £ s with minimum x-coordinate such that 
r is visible from p' . It follows that p' and r are collinear with a vertex q of some segment 
s' G S and that q is on the segment p'r (see Figure [5]a). Then q is an endpoint of an active 
edge e of V^{p') with start(e) to the left of p'. Since every active edge eventually adds a 
triangle to Cg{s), there is some pi to the right of p' that adds a triangle Ag to Cg{s) that 
contains r (see Figure [5}b). Since this is true for every point r G V^{s), we conclude that 
UC^(s) D V^{s), and hence UC^(s) = V^{s). □ 

Lemma 2. Let S be a set of n disjoint line segments. The covering Csis) can be computed 
in 

1. 0{ms) time if we are given EVG{S) or 
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2. Oe(n + (msra)^/^) otherwise. 

Proof. Part 1 of the lemma is clear. The algorithm for constructing Cs{s) processes the 
edges of EVG{S) incident on s in the order in which they appear. These edges can be 
easily extracted from EVG{S) in the order in which they appear and processing each edge 
takes 0(1) time. 

Part 2 of the lemma requires some use of a geometric range searching structure for 
answering ray-sweeping queries. Let q and q' be two points that are visible, with q' on some 
segment s" £ S. A ray-sweeping query asks to determine the first endpoint of a segment in 
S that is intersected by qq' as the point q' moves towards the left endpoint of s" . 

A ray-sweeping query is an optimization problem. It's corresponding decision prob- 
lem is a triangle interference query, which asks to determine if a query triangle A with 
vertices q, q' and q" S s" intersects any segment of 5. Because q and q' are visible and q' 
and q" are both on s" , it is not hard to see that if A does intersect some segment in S, 
then A contains an endpoint of a segment in S. That is, a triangle interference query can 
be solved using a triangular range searching structure built on the endpoints of segments 
in S. 

Triangular range searching is a well studied problem, and a number of solutions exist 
that, for any k with n < k < n^, give Oe{k) space structures with Of:{n/y/k) query time 
[U Section 4]. Using one of these structures and applying Chan's randomized optimization 
technique [HI Theorem 3.2] yields a data structure for ray-sweeping queries with the same 
preprocessing, space, and query time bounds. 

To construct Cs{s) we use essentially the same sweeping algorithm used to define 
Cs{s) except that ray-sweeping queries are used to compute the algorithm's events on the 
fly. The algorithm uses a priority queue Q to order and process these event points in left 
to right order. To initialize the algorithm, we construct the visibility polygon Vs{p = pi) 
in O(relogn) time using a radial sweep [UlIS]. Next, each active edge of Vs{p) is identified 
and processed. 

Anytime (during initialization or later) that an active edge e = qq' is created, the 
algorithm performs a ray-sweeping query with the segment qq' and a ray sweeping query 
with the edge q'p (see Figure [6]) . The results of these two queries determine an event point 
p' £ s to the right of p at which time the edge e contributes a triangle Ag to s. This event 
point p' is enqueued in Q. It is not hard to verify that this algorithm computes the same 
set of triangles Cs{s) as the original algorithm and that the number of ray-sweeping queries 
performed is O(m-s) (at most two queries are performed for each triangle added to Cs{s)). 

Therefore, the algorithm can be implemented to run in Oe{k + mgn/^/k) for any 
n <k <n^ . Given the value of rUg in advance, setting k = (mgn)'^/'^ would yield the stated 
time bound. However, even without knowing in advance we can begin by estimating 
the value of nig as m'^ = 2 and doubling our estimate (and rebuilding the ray-sweeping 
structure) if we discover that > m'g. This doubling strategy yields the overall time 
bound of Oe(n + (msn)^/^), as required. □ 

Next we show that, in a global sense, the number of triangles containing a point 
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Figure 6: The four cases that can occur when using a ray-sweeping query to determine the 
point at which active edge e = gg' contributes a triangle to Csis). 




Figure 7: The three cases that occur when processing an edge of EVG{S) incident on pi. 
Here, start(e) = pj. 

Lemma 3. Let C{S) = Uses ^s{s) and letp he any point in that is not on the boundary 
of any triangle in C{S). Ifrup is the number of segments in S (partially) visible from p and 
m'p is the number of triangles in C{S) that contain p, then rup < m'p < 2mp. 

Proof. Let Cp C C{S) be the set of triangles in C{S) that contain p, and let Sp C S 
be the set of segments in S that are (partially) visible from p. Our goal is to show that 
\Sp\ < \Cp\ < 2\Sp\. The lower bound on m'p = \Cp\ is trivial: For every segment s G Sp, 
Vs{s) contains p, so, by Lemma [ll Cs{s) contributes at least one triangle to Cp. 
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(a) /(A) = s (b) /(A) = s' 



Figure 8: The mapping / takes A onto (a) s and (b) s'. 

To prove the upper bound, we describe a mapping f : Cp ^ Sp that is 2-to-one; for 
every s G Sp, there exists at most two triangles A £ Cp such that /(A) = s. The existence 
of / then proves the upper bound. 

Let A E Cp be some triangle that contains p and suppose that A G Cs{s) for some 
s £ S that is, without loss of generality, below p. If A is incident on s (Figure [8ja), then 
A was added to Cs{s) as part of Vs{p) where p was the left endpoint of s. In this case, we 
set /(A) = s. Otherwise, A was created when sweeping s with p and some active edge e 
of Vs{p) generated A (Figure [8|b). The vertex g of A that is closest to s is incident on a 
segment s' G S. In this case /(A) = s' . 

We now argue that / is 2-to-one. Let s £ S he some segment and suppose, without 
loss of generality, that p is above s. Consider a triangle A G /^^(s) and observe that, by 
the definition of /, A has a vertex that is an endpoint of s. 

Note that there is at most one triangle in Cs{s)r\Cp that maps to s, and this triangle 
exists precisely if p is visible from the left endpoint of s. All that remains to show is that 
there is at most one additional segment s' G S, s' ^ s such that Cs{s') contains a triangle 
A with /(A) = s. 

Let A be such a triangle and suppose that A is incident to the endpoint q of s. 
Refer to Figure [9j The triangle A was generated by an active edge when processing s' . In 
particular, there is a subsegment pjpi C s' such that an active edge e of V^{p) sweeps over 
A when p travels from pj to pi. (Note, pj = start (e).) This implies that pi and pj are below 
s. Since p travels from left to right along e, this implies that q is the right endpoint of s 
because, otherwise, e would not be an active edge of Vg~{p). 

Thus far, we have established that at most one triangle in f~^{s) is incident to the 
left endpoint of s. To see that at most one triangle (A, discussed above) is incident to the 
right endpoint of s, suppose by way of contradiction that there are two such triangles A 
and A' with A G C's(s') and A' G Cs{s"). Consider the line £ through p and q. Observe 
that i intersects both s and s' , in two points r and r' , respectively. But this is not possible 
since then one of r or r' does not see the endpoint q. □ 

Remark: The condition, in Lemmajsj that p is not on the boundary of any triangle in Cs{s) 
is unnecessary if we take a little extra care. In particular, the mapping / actually maps 
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Figure 9: At most one triangle in / ^(s) is incident to the right endpoint of S. 

triangles to the endpoints of segments. The set Cs{v) of triangles mapped to a particular 
endpoint v all have w as a vertex and no two triangles in Cs{v) share an interior point. 
This means that we can define each triangle in Cs{v) to either include or exclude some of 
its edges or vertices so that the triangles are disjoint but their union remains unchanged. 
This yields a set of (partially open) triangles C'{S) for which Lemma [s] holds for any point 

3 Applications 

In this section, we consider applications of Lemma[T]and Lemma|3]to some visibility testing 
and counting problems. These applications rely on data structures for triangle inclusion 
counting: Given a set T of triangles, we want to preprocess T into a data structure for 
counting the number of triangles in T that contain a query point p. 

The tools needed to perform these queries are well-known, but finding the relevant 
structures and techniques, and applying them correctly, can take some time. Therefore, we 
review the data structure here and point out the relevant references. 

Let A be a triangle. Then A is the intersection of at most 4 halfplanes bounded 
by four lines /i(A) = [ui,U2-,di^d2) where ui and U2 bound A from below and di and d2 
bound A from above. Given a triangle A we have either ui = U2 or di = d2- By the 
standard duality mapping [71, Section 8.2] the four lines in /i(A) map to four points u\, 
d\ and d^- A point p G maps to a line p* in the dual plane. The point p is contained 
in A if and only if the line p* is above (or on) u\ and U2 and below (or on) d\ and d^- Let 
h*{^) = {ul,ul,d\,dl). 

A triangle inclusion counting structure for T stores the 8-dimensional point-set 

h*{T) = {/i*(A) : A G T} . 
Given a query point p, we want to count the number of points (a, 6, c, d) G h*(T) that satisfy 
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the four requirements: 

1. a is above p* , and 

2. 6 is above p* , and 

3. c is below p* , and 

4. d is below p* . 

Counting the number of points in h*{T) that satisfy any one of these requirements is a 
halfplane range counting problem. Data structures for halfplane range counting are plen- 
tiful, and there are several data structures known that use 0^{k) space and have query 
time Oe{n/y/k) [U Section 4]. Several of these structures (for example, Matousek's efficient 
partition trees ) are hierarchical structures that are efficient and r-convergent (see Agar- 
wal and Erickson Section 5] for definitions of hierarchical, efficient, and r-convergent). 
This implies [U Theorem 10] that there exists a 4-layer structure that uses Oe{k) space and 
preprocessing time and that, in time Oe{n/y/k), can count the number of elements in h*{T) 
that satisfy the constraints 1-4 for any query point p* . Translating this back into primal 
space we obtain the data structure we need: 

Theorem 1 ([UIIl]). Let T he a set ofn triangles. For any k with n < k < n?, there exists 
a data structure of size Oe{k) that can be constructed in time Oe{k) and that can count the 
number of triangles containing a query point p in Oe{n/^/k) time. 

3.1 Visibility Testing 

Our first application follows immediately by storing the triangles of Lemma [T] in the data 
structure of Theorem [T] This yields our first result: 

Theorem 2. Let S be a set of n disjoint line segments and let s £ S be a special seg- 
ment. For any nis < k < ml, there exists a data structure of size Oe{k) that can test, in 
Oe{ms/ \/~k) time, if any query point p is contained in Vs{s). The data structure can he 
constructed in 

1. Oe{k) time if we are given EVG{S) or 

2. 0^{n + (msTi)^/^ + k) time otherwise. 

Next we argue that, barring a breakthrough on Hopcroft's Problem j8j. Theorem [2] 
is near-optimal. Hopcroft's Problem takes as input a set L of n lines and a set P of n points 
and asks if any point in P is contained in any line in L. Currently, the most efficient methods 
of solving Hopcroft's Problem have running times in i7(n^/^). Furthermore, Q{n^^^) is a 
lower bound for Hopcroft's Problem in a restricted model of computation that can model 
all known algorithms for the problem [8]. 

Given the set L, we can compute the leftmost intersection point between any pair of 
lines by sorting the lines by slope and checking the intersection points between consecutive 
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(a) (b) 

Figure 10: A set L of lines (a) and a set of 3n + 0(1) segments where testing if a point is 
in Vs{sq) helps to determine if the point is contained in any line of L. 



pairs of lines. Assume, without loss of generality, that this leftmost intersection point has x- 
coordinate equal to 0. Using infinitesimal gaps between segments]^ we can easily construct 
a set of 3n + 0(1) segments so, ... , •S3n+o(i) such that a query point p whose x-coordinate is 
greater than is visible from sq if and only if p lies on one of the lines in L (see Figure 10). 
For a query point p with x-coordinate smaller than we can test if p is contained in any 
line of L in O(logn) time by storing the lines of L sorted by slope and using binary search. 

Therefore, by the above discussion, setting k = n^/^ and using Theorem [2] we can 
use this data structure to solve Hopcroft's Problem in Ot{n^^^) time. Furthermore, the 
existence of a data structure for testing if Vs{s) contains a query point p that could be 
constructed in o(n^/^) time and whose query time is o(n^/^) would give a o{n'^/^) time 
algorithm for Hopcroft's Problem. 



3.2 Visibility Counting - Relative Approximation 

Next we consider Fischer et a/.'s problem of approximate visibility counting |101 lllj . We 
want to preprocess the segments in S, so that for any query point p we can quickly approx- 
imate the number of segments in S that is visible from p. 

We begin with an easy corollary obtained by computing C{S) using Lemma [2] and 
putting all the triangles of C{S) into the data structure of Theorem [l| The resulting 
structure guarantees a relative approximation of the visibility count for all values of nip: 

Corollary 1. Let S he a set of n disjoint line segments whose visibility graph has m edges, 
and let < a < 1 be real valued parameters. There exists a data structure D that can 
approximate the number of segments of S visible from any query point p such that: 

1. D has size Oe(mi+") = Oe(n2(i+°)), 

2. D can be constructed in time Oe{m^^°') = Oe{n'^^^^°'^), 

^The use of infinitesimals in lower bounds is justified by Erickson's results [9]. 
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3. D can perform a query in Oe(m^^/^^^^ = Oe{n^ °) time, and 

4. when querying D with a point p that sees rup points of S, D returns a value m'p that 
satisfies nip < m'p < 2mp. 

3.3 Visibility Counting - Absolute Approximation 

Although Corollary [1] offers a good approximation guarantee, the space requirement is too 
large. In the worst case, when m = Q(n'^), a data structure of size (^{n"^) is required in 
order to achieve a sublinear query time. 

Fischer et al. [101 [TT] argue that, for the computer graphics application they consider, 
an absolute approximation is sufficient. In their application, there is a function f{n) such 
that, for mp <^ fin) it is more efficient to run a visibility culling algorithm before rendering 
the view from p but for f{n) <^ mp it is preferable to simply send all elements of S to the 
graphics hardware for rendering. For mp ~ fin) neither strategy has a clear advantage. If 
we define a<^b as a<b — 5n then we see that an algorithm that can approximate mp with 
an additive error of at most 5n is sufficient for this application. 

We present two different data structures that offer this kind of approximation guar- 
antee. These two structures offer different tradeoffs in terms of accuracy and query time. 

3.3.1 Solution 1: Sampling from S 

The data structure of Theorem [2] combined with a careful random sampling of the elements 
of S provides our first solution. We create a Bernoulli sample S" ^ S* by choosing each 
element of S independently with probability (c log n) /n, where c > 1 is a parameter of the 
data structure that controls the accuracy of the approximation. For each sample s G S" , we 
construct the data structure of Theorem [2] with the value k = mg"*"" for some parameter < 
a < 1 that controls the space /query-time tradeoff. If, during the construction of this data 
structure, it turns out that ms > 4:cn, then discard s from S" . Notice that this algorithm is 
effectively drawing a Bernoulli sample from the set S' = {s £ S : mg < 4cn} and that, since 
2m = < 4n^, there are at most n/c elements in S that are not in S' . Suppose p 

is visible from mp elements of S and m'p elements of S' . Then, by the above discussion, we 
have mp — n/c < m'p < mp. Let m'p = {n/{c\ogn)) - {{s £ S" : p £ Vs{s)}\. The quantity m'p 
is an unbiased estimator of m'p and, using Chernoff's bounds (see Appendix |A]) , we readily 
establish that 

Pr{|m^' - m'p\ > 6n} < ^-^(^'^/S) < n-^iS^^n/m^) 

for any 5 > 0. Combining this with the previous equation gives 

Pr{|m;,' - mp\ > 6n + n/c} < n-^('5'^"/™p). 

This establishes the accuracy of the data structure. What remains is to analyze the query 
time, space, and construction time. 

Query time. A query computes m" by performing a query in each of the data structures 
built on the elements of S" . The expected contribution of an element s £ S' to the query 
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time is therefore 

(^) • 

since it contributes Oe{rns/Vk) = 0^{mi^^'^^^^ to the query time if it is chosen to take 
part in 5"' and it contributes nothing otherwise. Summing this over all s, we get a total 
expected query time of at most 

X a (m(i/2)(i-")) = 0,(c(m/n)(^/2)(i-")) 
seS' 

where the last step follows from the fact that f{x) = x^^^'^^^^~°'^ is a concave function and 
that X^sgg/ mg = 0{m). 

Space. Arguing as above, the expected amount of space that an element s £ S' contributes 
to this data structure is 

(^^).o,w«). 

Therefore, the total expected amount of space used by the structure is 

f — ) • E = 0,((cm/n)(cn)") (1) 

where the last step follows by maximizing the sum YlseS' '"^s"''" using the facts that X]sg5' 
0{m) and that any individual s €z S' has rrig < 4cn. 

Preprocessing time. The preprocessing phase requires computing Cs{s) for each sample 
element s G S" and constructing a layered partition tree for the elements of Cs{s). Con- 
structing the partition tree takes Oe(mJ+") time, so, as above, the total expected cost of 
constructing the partition trees for all elements in S" is Oei{cm/n){cn)'^). 

Computing Cs{s), using Lemmajijtakes Oe((msn)^/^) time. Since /(x) = (xn)^/^ is 
a concave function, the total expected time to compute Cs{s) for each s G S" is Oeicm^^^). 

Theorem 3. Let S be a set of n disjoint line segments whose visibility graph has m edges 
and let c > 1 and < a < 1 6e real valued parameters. There exists a data structure D that 
can approximate the number of segments of S visible from any query point p such that: 

1. D has expected size Oe{{cm/n){cn)'^) = Oe{{cn)^~^°'), 

2. D can be constructed in Oe{cm?/^ + {cm/n){cn)"') = Oe{cn^/^ + (cn)^+") expected 
time, 

3. D can perform a query in Oe{c{m/ n)^^^'^^^^^"'^) = Oe{cn^^^'^^^^^"^) expected time, and 

4- for any 6 > 0, when querying D with a query for a point p that sees nip points of S, 
D returns a value m'p that satisfies rrip — nj c — bn < m'p < rrip + 5n with probability 

at least i _ n-f^('^'^"/'"p). 
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Figure 11: The two cases in which the ordered pair {u,v) generates the triangle A. 



Example: For any constant 6 there exists a c = c{5) such that taking a = 4/3 gives a data 
structure of size Oe{{m/n)n^/^) = 0^{n'^^^) with query time Oe((m/n)^/'^) = 0^{n^^^) and 
the structure approximates nip for any p with an absolute error of at most 5n w.h.p.. 

As this example shows, the structure of Theorem [3] is quite efficient for constant 
values of 6. Unfortunately, the theorem becomes weaker when using subconstant values 
of 6. This is because, to obtain meaningful error bounds, we require c = Cl{l/6) and the 
running time of the query algorithm grows linearly with c. 

3.3.2 Solution 2: Sampling C{S) 

Next we consider a different data structure that is also based on random sampling. Rather 
than sample segments of S, we instead sample triangles of C{S) and use Lemma [s] to bound 
the quality of the approximation. This results in a more efficient space/accuracy tradeoff 
than that provided by Theorem [3} The cost of this savings in space is that we obtain a 
relative approximation bound when rrip is large and an absolute approximation bound when 
rUp is small. For the application proposed by Fischer et al. [TT] this is an acceptable 
approximation bound. 

Consider the set C{S) of triangles described in Lemmajsj For any point p S M^, the 
number, m'p, of triangles in C{S) that contain p is a 2-approximation of the number, rUp, 
of segments in S that are visible from p. In particular 

rUp < m'p < 2mp. . (2) 

Our strategy is to approximate m'p by sampling elements of C{S). The easiest 
way to proceed would be to select a Bernoulli sample by sampling each element of C{S) 
independently with probability (c log n)/n. This would require enumerating the elements of 
C{S), of which there are 0(m), yielding a construction time that is O(n^) in the worst case. 
Instead, we use a different sampling strategy based on the rejection method that avoids 
computing C{S). 

Sampling C{S). Our goal is to obtain a random multiset C"{S) C C{S) of size roughly 
c{m/n) logn. To achieve this, we repeat the following procedure 4cnlogn times: We select 
two points u and v at random, with replacement, from the 2n endpoints of S. Note that 
there are 4n^ ways of doing this. Next, using 0(1) ray sweeping queries, we determine if q 
is a vertex of some triangle A G C (S) that has an edge collinear with wb and that lies to 
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the left of ut (see Figure [TTj) . Note that, for any A S C{S), there is exactly one pair {u,v) 
for which this is truej^ Therefore, if this test is affirmative then A is an element drawn 
uniformly at random from C{S) and we add it to our sample C"{S). The probability that 
we increase the size of C"{S) this way is ms/i^n^), where ms = \C{S)\ = 0{m). 

Space, preprocessing time, and query time. To compute C"{S) efficiently we use the 
ray-sweeping data structure described in the proof of Lemma[2j Each sampling step requires 
0(1) ray-sweeping queries, which can be done in Oe(n/\/Z) time after Oe{i) preprocessing. 
Thus, the expected time required to build the ray-sweeping data structure and perform 
4cnlogn sampling steps is Oe{i + {cn log n){n/Vi) = 0^(1 + cn'^ /Vi) = Oe(c^/^n^/^) for 
£ = c2/3n4/3. 

Each sampling step adds an element to C"{S) with probability m5'/(4n^). So, the 
number of samples in C"{S) is a binomial random variable with parameters 4cnlogn and 
ms/{4:n?) and the expected size of C"{S) is therefore c{ms/n)logn. Using Chernoff's 
Bounds, we find that the probability that the size of C"{S) exceeds ac{ms/n)logn is at 
most for any a > 1. This concentration result ensures that when building the data 

structure of Theorem [T]on the elements of C"{S) the expected size, preprocessing time, and 
query time of the resulting structure are Oe{k), Oe{k), and Oe{{m/n)/y/k), respectively, for 
any m/n < k < {rajn)^ . 

To summarize, for any k with m/n < k < {m/n)^, the above sampling procedure 
runs in O^^c^/'^n^/^ + k) expected time and produces a data structure of Oe{k) expected size, 
that can answer queries in Oe{{m/n) /^/k) expected time. All that remains is to calibrate 
and check the accuracy of the results provided by the data structure. 

Estimating m'^. Recall that our goal is to estimate m^, the number of elements of C{S) 
that contain the query point p, as m'^ is a 2-approximation to the number of segments of S 
visible from p. Let 

ml = (n/(clogn)) • |{A G C" {S) : p G A}|. . 

(Note that computing m'p does not require knowing the value ms = \ C (S) | .) Each step of the 
sampling procedure finds an element A G C{S) such that p £ A with probability exactly 
m'p/ (An?). Since the sampling procedure runs for 4cnlogn steps, this implies that the 
number of triangles in C"{S) that contain p is a binomial random variable with parameters 
4cnlogn and mp/(4n^). Therefore, 

E[mp] = (n/(clogn))(4cnlogn)(mp/4n^) = . 

That is, m'p is an unbiased estimator of m'p. Furthermore, applying Chernoff's bounds to 
the underlying binomial random variable (see Appendix [B]) , we find that 

Pr{|m;' - m'p\ > 6n} < ^-^^(5'^".^) 

for any 5 > 0. Combining this with ^ we obtain 

Pr{mp -5n< m'^ < 2mp + 5n} > 1 - ^-^('^'^'"f) . 

''The pair {u, v) generates the triangle A precisely if u is below v and v is the right endpoint of its segment 
or u is above v and v is the left endpoint of its segment. 
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This establishes the accuracy of the data structure and completes the proof of our last 
theorem: 

Theorem 4. Let S be a set of n disjoint line segments whose visibility graph has m edges 
and let c > 1 and < a < 1 6e real valued parameters. There exists a data structure D that 
can approximate the number of segments of S visible from any query point p such that: 

1. D has expected size Oe{{cm,/nY^"^) = Oe((cn)^+")), 

2. D can be constructed in time Oe{<?^^n^^^ + (cn)^"*"°), 

3. D can perform a query in 0^{{cm/n)^^/'^^^^^°'^) = Oe{{cn)^^^'^^^^~"^) time, and 

4- for any 6 > 0, when querying D with a point p that sees nip points of S, D returns 
a value m'p that satisfies nip — 6n < m'p < 2mp + 6n with probability at least 1 — 

^—fl(S^cn/rnp) 

Example. Taking c = dn^^^, for a large constant d, and a = 0, we get a data structure 
of size 0,(n4/3) that can be constructed in time Oe(n^^/^) and that can, in Oe(n^/'^) time, 
effectively distinguish between viewpoints p where rup <C n^/^ and viewpoints p where 
nip v?!'^ . 

4 Summary and Conclusions 

Many open questions remain. The data structure for testing if a point is in 1/5(5) for 
a segment S (Theorem [2| is near-optimal, at least assuming an 0(n^/^) lower-bound for 
Hopcroft's Problem. However, it is difficult to say if the data structures for approximate 
visibility counting are close to optimal. Our solutions reduce visibility counting to the 
problem of computing the depth of a query point in an arrangement of 0((m/n) log n) 
(Theorem [3] and Theorem |4]) or Oira) (Corollary [T]) triangles. In both cases, it would be 
sufficient to give a relative approximation for the depth of the query point. Unfortunately, 
without some additional assumptions (such as fatness) about the triangles, there is currently 
no good solution to this problem. 

The results in the current paper consider the problem of planar visibility counting, 
where S is a set of disjoint line segments in M?. Of course, modern virtual environments are 
often 3-dimensional. Many of these environments are just barely 3-dimensional in the sense 
that they consist of a constant number of 2-dimensional layers that can be handled using 
the data structures presented in the current paper. However, ultimately we would like to 
develop data structures that store a set S of disjoint triangles in M'^ and can approximately 
count the number of elements of S (at least partly) visible from a query point p E M^. 
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Appendix 

A Accuracy Bound for Theorem |3] 

In this appendix, we derive the error bound of the data structure of Theorem |3j To do this, 
we will use a version of Chernoff Bounds for binomial random variables [21 Appendix A.l] 
which states that, for a binomial random variable B with mean /i, 

Vi{\B - ii\> r/i} < exp(-0(rV)) • (3) 

for any r > 0. 

Let B be the number of samples in S" visible from p, let x = m'p be the number of 
segments in S' visible from p, and let t = (clogn)/n. Then B is a binomial(a;, random 
variable with expectation fi = xt. We have that \mp — m'p\ > 5n if and only ii\B — ^\ > t6n. 
Taking r = t5n/ ji and applying Equation ([s]) we obtain 



Pr{|mp — m!p\ > dn} = Pr{\B — fi\ > t5n} 

< exp 

= exp(-fi((Mn)VM)) 

= ex.p{-n{{6n)H/x)) 

= exp(— r2((5^cnlogn)/x)) 

= n~ 



-n{6^cn/x)) 



as required. 
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B Accuracy Bound for Theorem |4] 



Let B be the number of sample triangles in C"{S) that contain p, let x = m'p be the number 
of triangles in C{S) that contain p, and let t = 2;/(4n^). Then B is a binomial(4cn log n, t) 
random variable with expected value ^ = 4tclogn = (cx log n)/n. 

We have that \mp — m'p\ > 5n if and only if /x| > ^clogn. Takingr = (5clogn)//x 
and applying Equation ([s]) we obtain 

Pr{|mp — m'p\ > 6n} = Pr{|i? — ^| > 6clogn} 

= exp(-17((((5clogn)//i)V)) 
= exp(— r2((5clogn)^//i)) 
= exp(— r2(5^(cnlogn)/x)) 

^ ^-n(5^cn/x)) 
^ ^-Q(52cn/m;)) 

as required. 
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